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Summary 


The  Fir it  Technical  Report  contained  a  detailed  simulation  study  about  the 
possible  accuracy  of  (utAonomic  positioning  in  the  field.  Namely  we  refer 
to  Chapter  Three,  especially  to  FiguA.es  2-29.  -The  simulation  study  was 
based  on  the  assumption  that  the  objects  of  field  observations,  the  stars, 

are  perfectly  known,  including  also  the  various  astronomic  reductions  such 

c 

as  proper  motion,  aberration  Darallax  etc.  Here  we,  therefore,  allow 
some  uncertainty  of  the  star  position:  We  present  a  detailed  analysis  of 
how  to  proceed  once  only  dispersive  prior  information  about  star  positions 
is  available.  The  nonlinear  condition  equations  are  linearized  by  a  two- 
step  Taylor  series  expansion  accounting  for  stochastic  information  a,  6  of 
right  ascension  and  declination,  the  parameters  of  star  positions.  A  nu¬ 
merical  analysis  will  be  presented  in  the  Final  Report.  .■ 


f 
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0.  Introduction 


Within  the  First  Technical  Repent  (March  1987)  we  have  characterized  astro¬ 
nomic  positioning  as  the  determination  of  the  direction  parameters  of  the 
gravity  vector,  namely  of  astronomical  longitude  A  and  as tronomical  lati¬ 
tude  <l> ,  relative  to  an  earth-fixed  reference  frame.  According  to  the  ob¬ 
servational  techniques  cf  stellar  objects  by  means  of  a  theodolite,  an 
additional  quantity,  the  orientation  unknown  £  in  the  horizontal  plane 
has  to  be  determined.  Within  a  GauH-Helmert  model  of  condition  equations 
with  unknowns  we  simulated  various  observational  configurations  in  order 
to  get  an  accuracy  estimate  of  the  field  operations  by  which  astronomical 
longitude  A,  latitude  $  and  the  orientation  unknown  z  are  determined. 

Namely  we  refer  to  Chapter  Three  of  the  Tint  Technical  Report  (March  1987), 
especially  to  Figure-i  2-29.  The  accuracy  of  the  observations  had  been  as¬ 
sumed  as  follows:  The  r.m.s.  value  for  the  horizontal  direction  was  sta¬ 
ted  as  o-j-  =  1",  for  the  vertical  direction  og  =  1"  and  for  time  measure¬ 
ment  oQ  =  0.1  sec . 

O 

The  simulation  study  was  performed  under  the  restriction  that  the  respective 
star  positions  are  given  perfectly.  In  reality,  the  star  positions,  namely 
right  ascension  a  and  declination  6,  have  to  be  characterized  by  a  variance- 
covariance  matrix.  Its  influence  on  astronomic  positioning  is  analyzed  here¬ 
in.  We  derive  observational  equations  in  which  the  uncertainty  of  the  star 
position  is  represented.  Only  through  such  a  general  model  can  the  accuracy 
of  astronomic  positioning  in  the  iield  be  real i sti call ly  estimated. 
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1 .  Relations  between  the  reference  frames  in  geodetic  astronomy 

First  of  all  in  this  chapter  the  fundamental  relations  in  geodetic  astronomy 
between  observations,  unknowns  and  given  coordinates  shall  be  derived  which 
will  be  needed  further. 

1.1  The  systematical  structure  of  the  reference  frames 

The  reference  frames  used  in  geodetic  astronomy  may  be  arranged  on  different 
levels  which  are  numbered  in  turn  or  indicated  by  symbols:  0  corresponds  to 
1  to  *,  2  to  •,  3  to  o.  One  fundamental  vector  V1  belongs  to  every  level  i. 
In  details  this  is  as  follows: 


V  =  V'  =  z 


v  =  v*  =  -r 


v  =  v  =  n 


v3  =  v°  =  * 


the  position  vector  from  the  point  of  observation  to 
the  target  object  (terrestrial  or  celestial),  which  is 
generally  a  star; 

the  negative  gravity  vector; 

the  earth  rotation  vector  (it  has  the  direction  of  the 
axis  of  the  earth,  points  to  the  North  Pole  and  has  the 
value  of  the  earth  rotation  rate); 

the  ecliptic  normal  vector  (it  points  to  the  northern 
pole  of  the  ecliptic). 


An  orthonormal  reference  frame  E1  belongs  to  every  level  with  its  base  vectors 


as  follows: 


=  norm  V1 


norm  (V1*-1  x  V1) 


E1  x  E1 

~2  x  h 


(1-1) 


(1-2) 


(1-3) 


Here  "norm"  denotes  the  abbreviation  for  the  normalization  of  a  vector,  and  "x1 
the  vector  product. 


i 


New  reference  frames  are  at  the  lower  end  the  observational  frame  E'  of  the  level 


-  4  - 


"0",  whose  third  base  vector  is  located  in  the  direction  of  the  observation  and 
which  is  unique  since  it  is  not  a  reference  frame  in  the  literal  sense,  because 
there  are  no  vectors  described  with  regard  to  this  frame,  and  at  the  upper  end 
the  ecliptic  frame  E°,  which  has  hardly  any  practical  importance. 

In  addition  to  the  systematical  E-triads,  F-triads  also  appear  on  each  level. 
These  systems  have  the  common  third  base  vector  with  the  appertaining  E-frame, 
nevertheless  the  direction  of  the  first  and  second  base  vector  does  not  follow 
from  the  systematic  structure  of  the  fundamental  vectors  V1,  V1+1,  but  from  a 
more  or  less  arbitrary  definition. 

A  new  F-triad  is  the  theodolite  frame  F*,  whose  first  base  vector  F  lies  in 

2r  a: 

the  direction  "zero"  of  the  azimuth  circle  of  a  theodolite  which  is  set  up  in 
the  astronomical  horizon.  The  longitudinal  angle  (see  below)  of  an  observed 
direction  in  the  local  horizon  frame  is  the  horizontal  direction  T  and  is  re¬ 
corded  systematically  in  the  clockwise  direction,  but  conventionally  counter¬ 
clockwise,  T  =  -T  .  The  latitude  angle  (see  below)  is  the  vertical  direction 

S  C 

as  in  the  horizon  frame  E*. 

The  transformation  from  a  frame  E1  to  the  appertaining  frame  F1  is  always  a 
counter-clockwise  rotation  round  the  common  third  axis  with  the  orientation 
angle  H1, 

F1  =  R3(Hi)Ei  .  (1-4) 


is  the  rotation  matrix,  which  describes  a  rotation  of  a  frame  round  the 
third  axis.  It  is 


R3(y) 


cosy  Siny 

-siny  cosy 

0  0 


0 

0 

1 


(1-5) 


Corresponding  to  eqn.  (1-5)  the  rotation  matrices  for  the  rotations  round  the 
first  and  second  axis  read 


w-a. 
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Rj(a) 


cosa  sina 


0  -sina  cosa 


(1-6) 


B2(3) 


cosB 


sinB  0 


-sinB 


cosB 


(1-7) 


R  ,  R2  and  R3  are  also  called  elementary  rotations.  The  orientation  angles  H1 
(see  eqn.  (1-4)  )  are  in  detail: 


H  =  H*  =  I 


H  =  H  = 


H3  =  H° 


Gr ,  s 


the  orientation  unknown  of  the  theodolite  which  has 
been  set  up; 

the  Greenwich  sidereal  time; 


the  angle  between  the  line  of  intersection  of  the  eclip¬ 
tic  with  the  mean  galaxy  plane  and  the  direction  to  the 
vernal  equinox  ±90°. 


ft 


For  the  transformation  from  a  frame  E1+1  to  the  underlying  frame  E1  one  needs 
the  longitudinal  angle  x^+1  and  the  latitude  angle  <E>^+1  of  the  fundamental 


vector  V1  with  regard  to  the  frame  E1+1 


r =  -r2<9°o-*;+1)  wx-+1) 


-i+1 


(1-8) 


-E^xi+1 


i+1 


Re  is  the  special  case  of  a  rotation  matrix  of  Eulerian  type,  in  which  the  three 
elementary  rotations  are  connected  in  a  row  as  follows: 


first  rotation  round  the  third  axis,  R3(y  ) 


! 


second  rotation  round  the  new  second  axis  with  the  angle  (90°-6) ,R2(90°-B) 


thirH  rntaHnn 


ble.  1 :  The  E  and  F  reference  f rames 


In  the  matrix  R  the  second  rotation  round  the  second  axis  would  take  place 
-e 

with  the  angle  B,R2(B). 

These  longitudinal  and  latitude  angles  are  in  detail: 

v°  =  v'  =  A  the  azimuth  of  the  observational  direction; 

A1  A*  s 

*0  _  *i  _  d  the  vertical  direction 

1 

v1  =  y*  =  0  the  sidereal  time 

A2  A.  s 

O1  =  $*  =  o  the  astronomical  latitude. 

2 


For  the  transformation  from  a  frame  E1+  to  the  frame  F1  lying  diagonally  under¬ 
neath,  one  needs  additionally  the  orientation  angle  H1  (see  Fig.  1): 


F1  *  5E(xtl.*i+1.H1)| 


i+1 


(1-9 


For  the  transformation  from  a  frame  Fi+1  to  the  frame  E1  lying  diagonally  under¬ 
neath,  one  needs  the  longitudinal  angle  A*+1  and  the  latitude  angle  4>*+1  of  the 
fundamental  vector  V1  with  regard  to  the  frame  F1+1  : 


I1  -  -W.i  •*(..>! 


i+1 


(1-1 


The  latitude  angles  are  the  same  as  above,  the  longitudinal  angles  are  in  de¬ 
tail: 


the  horizontal  direction  of  the  observation  direction, 
systematically  measured  counter-clockwise,  conventionally 
in  the  clockwise  direction,  T  =  -T  ; 

s  c 


/r  =  a*  =  a 


the  astronomical  longitude; 


a:  =  A  =  90° 


3 

and 


,2 


=  4>  =  90°-t 


the  orthogonal  complement  to  the  inclination  of  the 


nrJVjvCTuvvTr^nr 
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level  5 


level  4,  mean  galaxy  frame 


Re(90°,e) 


level  3  (°),  ecliptical  frame 


-R3(6Gr,s> 


level  2  (°),  equatorial  frame  (fixed  in 
space,  fixed  in  the  earth) 


tL(6  ,$) 

-E  s 


Re(a/v) 


R3(i) 


level  1  ('),  horizontal  and  theodolite 

frame 


R  (A  ,B) 

“E  s 


Rf(t  ,b) 

E  s 


8,(0 


1  Commutative  diagram  with  reference  frames  in  geodetic  astronomy 
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For  the  transformation  from  a  frame  F1+1  to  the  underlying  frame  F1,  one 
needs  additionally  the  orientation  angle  H1: 


ZL  -  MAi+i’*i+i’Hi)f 


i+1 


(1-11 


The  transformation  in  the  opposite  direction  results  from  the  respective  trans¬ 
posed  of  the  rotation  matrix. 


1 .2  The  commutative  diagram  and  the  fundamental  equation  of  geodetic  astronomy 

In  the  rotations  introduced  up  to  this  point  the  right  ascension  a  and  the  de¬ 
clination  6  are  still  missing;  these  describe  the  observation  direction  E3, 
to  a  star  with  regard  to  the  space  fixed  equator  frame,  namely  in  the  same  way 
as  A  and  B  do  this  with  regard  to  the  horizontal  frame.  Indeed,  the  rotation 
FL(a,6)E*  does  not  lead  to  the  frame  E';  in  E"  =  FL(a,6)E  the  second  case  vec- 

—  L  ^  —  Cj  ~ 

tor  E2„  lies  in  the  equatorial  plane,  but  in  E'  the  vector  E2,  lies  in  the  hori¬ 
zontal  plane  (east).  Of  course  E'  and  EM  have  the  common  third  base  vector, 
but  they  differ  by  a  rotation  round  this  vector  with  an  angle  C  : 

E"  =  R3(C)E'  .  (1-12 

A  serial  connection  of  several  transformations  is  called  a  diagram  in  the  alge¬ 
bra.  If  the  transformations  which  have  to  be  reversibly  unequivocal,  form  a 
closed  circle,  this  is  called  a  commutative  diagram.  Then  one  is  able  to  express 
a  transformation  by  means  of  the  others.  Such  a  commutative  diagram  is  presen¬ 
ted  in  Fig.  1  by  the  lines  which  are  thickly  marked.  For  example  the  trans¬ 
formation  E*  -*■  E'  can  be  expressed  by 

E'  =  RE(T,B)R3(Z)E*  (1-13 

or 

E*  =  R3(C)RE(a,6)R3(0Gr)RE(A,'.‘')E*  . 

As  the  representation  of  the  triad  E1  with  regard  to  the  triad  E*  is  unequi¬ 
vocal  ,  it  follows  that: 


1 

i 

» 


(1-14 


RE(T,B)R3(Z)=R3(C)RE(a,6)R3(0Gr)RE(A,4.) 

This  equation  reads  with  the  right-hand  side  written  in  full  length 
RE(T,B)R3(l)=R3(-f;)R2(90o-6)R3(a)R3(-0Gr)R3(-A)R2(<I>-90o) 

and  can  again  be  reduced  to 
RF(T+E,B)=RTU,6-90°)R'I(er  +A-«,i) 

-b  -e  E  Gr 


These  are  the  desired  ^undarmntal  -*i elation*  between  the  parameters  appearing 
in  geodetic  astronomy. 
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2.  The  observation  equations  of  geodetic  astronomy 

The  fundamental  equation  (1-17)  consists  as  a  matrix  equation  of  nine  sepa¬ 
rate  equations,  of  which  only  three  are  independent  of  each  other  because  of 
the  property  of  orthonormal ity  of  the  rotation  matrices.  These  three  inde¬ 
pendent  equations  represent  concLction  equations  uxith  unknowns  for  every  star 
(if  T,B  and  0„  are  measured  at  one  instant). 

The  matrices  on  the  left  and  right-hand  side  of  equation  (1-17)  read  as  fol¬ 
lows  when  they  are  multiplied  respectively: 


cos(Z+T)sinB 


-sin(Z+T) 


cos(Z+T)cosB 


sin(Z+T)sinB 


cos(Z+T) 


sin(l+T)cosB 


Column  1 


_  .n$cos(0  +A-u)sin6cos£-sin<1sin(6_  +A-a}sim';+coS‘1)cos6cosf; 
Gr  Gr 

sinOcos(  +A-a)sindsinC+sin4sin(  6_  +A-a)cosi+cos1'COS^sinC 

Gr  Gr 

sintcos(C<  +A-a)cos6-cosvsini 
Gr 

Column  2: 


-sin(G^,  +A-a)sin6cosC-cos(0_  +A-a)sin£; 

Gr  g, 

-sin(0„  +A-a)sinosin£+cos(  £•  +A-u)cos-A 

Gr  Gr 

-sin(0^  +A-a)cos6 
Gr 

Column  3: 


cos<J>cos(CL  +A-ct)sin6cosC-cos:isin(i-  +A-  0  s i  n.-.-s  i  n :  cosZcosC 

Gr  Gr 

cosOcos(0_  +A-a)sinijSin£>C0S'tsin(C'  +A-a)coS‘>sinvCOS5sin£, 
Gr  Gr 


cost’Cos(  8_  +A-a)cos5+sinTsin6 

Gr 


WWHiJ  ■  UJJ 
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In  the  right-hand  matrix  the  elements  in  the  third  row  are  the  shortest  and  at 
the  same  time  the  only  ones  which  do  not  contain  the  angle  £.  Therefore,  it 
is  the  obvious  thing  to  do  to  select  two  equations  from  this  row  as  indepen¬ 
dent  equations.  As  a  third  equation  one  could  take  an  element  from  another 
row  of  the  matrices  whereby  the  angle  £  ,  in  which  one  is  not  actually  in¬ 
terested,  would  indeed  appear  as  an  additional  unknown.  So  one,  therefore, 
dispenses  with  such  an  equation  and  there  remain  only  two  independent  equa¬ 
tions  for  one  complete  observation  (T,B  and  b  )  with  the  three  unknowns  A, 

Gr 

4>,I: 


sinB  =  cos<tcos(0^  +A-a)cos<5+sin$sin{ 

Gr 

(2-1) 

sin(Z+T)cosB  =  -sin(0„  +A-a)cos6 

Gr 

(2-2) 

cos(I+T)cosB  =  sinicoste,,  +A-a)cos4-cosOsin6 

Gr 

(2-3) 

Equations  (2-1)  and  (2-2)  are  independent  of  each  other,  equation  (2-3)  is 
dependent  on  them  both.  It  will  be  used  later  only  for  the  determination  of 
approximate  values.  The  appearing  variables  be  summarized  once  more: 


A 

* 

a 

h  =  0  +  A  -  a 

Gr 

E 

T 

=->A  =  E  +  T 
B  =  90°  -  z 


astronomical  longitude 
astronomical  latitude 
right  ascension  of  the  star 
declination  of  the  star 
hour  angle 

orientation  unknown  of  the  instrument  (theodolite) 

horizontal  direction;  observed 

azimuth 

vertical  direction,  angle  between  horizon  and  star;  observed 
Greenwich  apparent  sidereal  time;  observed 
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2.1  Linearization  and  matrix  representation 

The  equations  (2-1)  and  (2-2)  are  linearized  by  Taylor  series  where  we  build 
in  pnA.cn.  stochastic  information  of  the  star  positions  a,  6. 

Observations: 

The  horizontal  direction  T,  the  vertical  direction  B  and  the  coordinate  clock 
time  t  are  obsenved.  In  order  to  derive  linearized  condition  equations  from 
(2-1)  and  (2-2)  horizontal  and  vertical  directions  are  decomposed  into 


T  =  t  +  At 


b  +  Ab 


The  coordinate  clock  time  x  has  to  be  related  to  the  Greenwich  apparent 
sidereal  time  angle  e^,  e.g.  by  the  series 


0Gr(t)  =  eGr(To)  +  V(to)(t  "  V  +  °2[(t  '  To}  1 
T  =  Tn  +  0r.^J  ~  0r..(O]  • 


(2-5) 


(2-6) 


0Gr  may  be  identified  with  the  earth  rotation  speed  n. 


Condition  equations: 


A nsatz : 


F [ E { Y } ,  Xp  X2]  =  0 


(2-7) 


The  nonlinear  vectorial  condition  equation  F[E{Y),  X^,  X^]  contains  the  ob¬ 
servation  vector  Y,  {stochastic  according  to  the  theory  of  measurements), 
the  fixed,  unknown  vector  X^,  but  with  dispensioe  stochastic  pnion  information 
being  available,  and  the  fixed  unknown  vector  with  any  stochastic  prior  in¬ 
formation  . 


1st  Taylor  series  ] 

At  a  fixed  point  y,  x^.  x^  of  approximation  we  are  linearizing  the  system 
of  nonlinear  condition  equations  (2-7). 


F[E{Y},  X 1 ,  X21  =  F(y,x1,xz)  +  F'ly.x^x^  [ E{ Y)  -  y] 


(2-8) 


+  Fx1(y’xl’x2,(Xrx1)  +  Fx2(y’xrx2)(X2-x2)  +  °2 


o0  indicates  the  neglected  terms  of  higher  order. 


I 
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2nd  Taylor  series 

At  a  stochastic  point  y,  of  approximation  we  are  linearizing  the  system 


of  nonlinear  condition  equations  (2-1). 


F(y,XpX2)  =  F{y,x1,x2)  +  F  '  (y ,x  1 , x2 )  (y-y ) 


+  F^  (y,x1,x2)(x1~x)  +  o2 


(2-8) 

(2-9) 


(2-9) 


F[E{Y},XrX2]  -  F(y,xrx2)  =  F‘ (y,x^  ,x2)  [  E  { Y }  -  y] 

+  ( y > x ^ , x2 ) ( X ^ -x ^ )  +  F^  (y , x^ , x2 ) ( X2~x2 ) 


(2-10) 


Now  we  have  been  led  to  a  system  of  tineaAized  condition  equations  where 
Y  -  y  =:  Ay  and  -  x^  =:  ax^  are  stochastic  collection  vectors  while 
X2  -  x2  =:  ax2  contains  the  fixed  or  non -stochastic  collection  oecton..  The 
next  step  is  to  write  down  the  obsenvationaJL  equation  for  y  and  the  pseudo- 
obsenvationa!  equation  for  x^,  the  vector  of  stochastic  information. 


E { Y }  ,  D { Y }  =  Qo  ,  rk  Q  =  n 

E{aj>  =  Xj  -  b  ,  DfXj}  ^  V^2  ,  rk  V2  £  m} 


(2-11) 


Here  the  pseudo- obsexvationaZ  equation  is  generated  by  prior  informative  data 
of  the  stcui  positions ,  parameterized  by  {a, 6},  namely  by 


E  {  .  !  }  = 

6  I 


,  D  {  '}  =  V.o  , 

6  ! 


x ^  =:  [a, 6] 

sinB  =  sin(b+Ab)  =  sinb  +  Ab  cosb 
cos*  =  cos(*  +6*)  =  cos*  -  6*  sin<l> 

O  O  0 

C0S(6Gr+A-a)  =  COs(09Gr+AQ-a+6eGr+6A-6a)  = 

=  cos^0Gr+Ao-a)  -  (60Gr+«A-«a)sin(oeGr+Ao-5) 
cos6  =  cos(6+66)  =  cos6  -  66  sin6 


(2-12) 


(2-13) 


(2-14) 


sin*  =  sin(*  +6*)  =  sin*  +  6*  cos* 


(2-1),  (2-2),  (2-14) 


.  ;  sinb  +  cosb  Ab  =  cos?  cos(  8r  +A  -a)cos6  +  sin*  sin6  -  (2-15) 

o  o  Gr  o  o 

-  cos<J>o  sin(oeGr+Ao-i)(606r+6A) 

I  +  [-sin*  cos(  e_  +A  -a)  +  cos*  si n6 ]  6* 

j  o  'o  Gr  o  '  o  J 

!  +  cos*  sin(  8r  +A  -a)  6a  + 

'  o  o  Gr  o 

I 

!  +  sin*Q(cos6)  66  ; 

sin(ZQ+t)  cosb  +  cos(lQ+t)  cosb(6E+At)  -  (2-16) 

-sin(z  +t)  sinb  Ab  = 

-sin(  0„  +A  -a)  cos6  -  cos(  Qr  +A  -a)  cos6(60r  +6A) 
o  Gr  o  o  Gr  o  '  Gr 

+cos(o0Gr+AQ-a)6a  +  si n ( 09gr+A0-a ) ( si n6 )  66 

In  extension  of  the  linearized  condition  equations  (2-4),  (2-5),  respectively 
being  introduced  in  the  Fin. 6 t  Technical  Rep onX,  page  13,  here  we  are  led  to 
the  additional  terms  proportional  6a,  66,  respectively.  They  account  for  the 
unceAtaintia  ci&  the  stoji  positions.  Once  we  assume  zero  F^.X^X^)  - 
F ( y » X ^ , X2 )  =  0,  a  result  which  can  be  achieved  by  a  proper  gauging  of  the 
approximate  values,  we  finally  introduce  the  linearized  equations 


-cosb  Ab  -  cos$  sin(  0„  +A  ~a)(60r  +6A)  + 

0  0  Gr  o  n  Gr 

+  [-sin*  cos(  9r  +A  -a)  +  cos*  sin61  6$  + 

L  0  '0  Gr  0  0  J 

+  cos<t>osi n(o0Qr+Ao_a) 6a  +  sin*Q(cos6)  66  =  0 


(2-17) 


-cos(ZQ+t)  cosb(6l+At)  +  sin(z  +t)  sinb  Ab  - 

-COS(o0Gr+Ao-a)  COS6 ( 60gr+6A)  + 

+cos(  0r  +A  -a)  66  +  sin(  0r  +A  -a)(sin6)  66 
'0  Gr  0  0  Gr  0 


(2-18) 


From  these  equations  we  obtain  the  proper  condition  equations  with  respect  to 
the  estimation  /  prediction  problem  when  we  introduce  the  vector  incon¬ 
sistency,  e^  and  e^  ,  respectively. 


(2-19) 


!  O 


•  A.b.  +  A-Ax„  +  BAy  =  A,e,  +  Be 
1  1  2  2  11  y 

A^,  A^  and  B  are  the  coefficient  matrices  for  all  observations  taken  to  various 
stars.  The  structure  of  the  coefficient  matrices  can  be  read  from  the  two 
equations  ( 2- 17 )  and  (2-18).  Note  that  not  only  A,  4>,  x  or  6A,  6$,  6X  are 
unknoi md,  but  as  well  a,  6  or  6a,  66,  the  corrections  to  star  positions. 

Their  uncertainty  will  play  an  essential  part  in  the  final  computation  of  the 
variance-covariance  matrix  of  A,  4>,  x! 


In  order  to  solve  the  linearized  mixed  model  condition  equations  (2-22)  we 
apply  the  g<zwAalize.d  method  o&  l<zoa,t- 


L{ey,ex  ,ax2,A} 
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(ey,e  ,ax2,a) 

(2-23) 


The  normal  equations  and  their  solution  will  be  given  in  the  Final  Report 
There  an  extension  of  the  results,  F iru>t  Technical  Rzpont,  pages  20-29, 
will  be  presented. 


